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ABSTRACT 

To date, there have been several detections of high-mass black hole binaries in both 
the Milky Way and other galaxies. For some of these, the spin parameter of the black 
hole has been estimated. As many of these systems are quite tight, a suggested origin 
of the spin is angular momentum imparted by the synchronous rotation of the black 
hole progenitor with its binary companion. Using Cygnus X-l, the best studied high- 
mass black hole binary, we investigate this possibility. We find that such an origin 
of the spin is not likely, and our results point rather to the spin being the result of 
processes during the collapse. 
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1 INTRODUCTION 

Observations have so far led to a large number of discovered 
high-mass X-ray binaries, both in the Milk y Way and othe r 
galaxies (over 100 in the Milky Way alone; iLiu et al.| [2006U 
However, only a handful of these are believed to harbor black 
holes. Two of these, Cygnus X-l and Cygnus X-3, are found 
in our own Galaxy. The most well-known and well-studied 
of the black hole binaries is Cygnus X-l. 

In parallel with the discovery of more sources, new tech- 
niques have given estimates of the spin parameter of the 
black hole. The two main techniq ues rely on spectr al fitting 
of the thermal component (e.g.. IShafee et ahll2006t) or rela- 
tivistically broadened iron line (e.g., lMilrerll2007l h However, 
these results are quite sensitive as spectral fitting is prone 
to a certain amount of degeneracy and model dependency. 
The technique is also dependent on the spectral state of the 
source, making it difficult to apply to certain sources which 
do not show this st a te, e.g . Cyg X-l. 

lAxelsson et al.l l|2005f ) presented an alternative tech- 
nique for determining the spin parameter in Cyg X-l. By 
studying the evolution of quasi-periodic oscillations (QPOs), 
they were able to give support to the idea that the oscilla- 
tions are connected to the relativistic precession frequencies 
predicted by general relativity in the strong gravitational 
field close to the black hole. Based on this identification, the 
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spin parameter was measured to be a* = 0.48 ±0.01, assum- 
ing the black hole mass to be 9 Mg. While still dependent 
on the theoretical interpretation, this method allows for a 
much more precise determination of the spin than the spec- 
tral modelling techniques. 

As several of the black hole systems discovered so far 
are tight syst ems, a possible source of sp i n may be tidal lock - 
ing (see, e.g.. iPodsiadlowski et aill2004l : iLevan et alj|2006h . 
Upon collapse of the black hole progenitor, the angular mo- 
mentum is preserved in the spin parameter of the black 
hole. In this paper, we will investigate whether such a sce- 
nario can explain the measured spin parameter of Cyg X-l. 
This system is a good candidate for such an investigation: 
it is bright and thereby we ll studied, and observations by 
iMirabel fc Rodriguesl (|2003h indicate that the mass loss in 
the formation of the black hole was low. 

In the rest of this section, we motivate our choice for 
the system parameters such as the mass of the compact ob- 
ject. Based on this, in the following section we investigate 
whether systems similar to Cyg X-l can acquire sufficient 
angular momentum via tidal locking to explain the current 
spin parameter. In Sect. [5] we discuss the evolution of the 
system in the past, including possible processes which may 
cause the spin measured today to differ from the natal one. 

1.1 The mass of the compact object 

Although one of the best studied black hole candidate 
sources, estimates of the current parameters of the Cyg X- 
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Parameter 



Value 



Reference 



Mass function (/) 0.251 M 

Black hole mass (Mco) 8-10 M Q 

Inclination (i) 28°-38° 

Separation (do) 32-42 R Q 

Donor mass (M 2 ,o) 7-21 M 

Period (P) 5.59982 d 

Spin parameter (a, ) 0.47-0.49 



Gies et al. (2003) 
See text for details. 
Gies fc Bolton (1986) 



From /, Mco and i. 
From /, Mco and i. 
Gies et al. (2003) 
Axelsson et al. (2005) 



Table 1. Table of parameter values used in this paper. The sepa- 
ration and donor mass ranges are calculated using the mass func- 
tion with the assumed ranges of black hole mass and inclination. 
Note that these ranges are not to be seen as confidence intervals 
in the traditional sense, but merely the values possible when M co 
and i are allowed to vary over our chosen ranges. 



1 system vary. Perhaps most well-determined is the mass 
function /, 



(Mco sini) 3 
(M co + M 2 ) 2 



(1) 



where M co and Mi are the masses of the compact object 
and compa nion star, respect ively, and i the inclination of 
the system. iGies et alJ l|2003r i found / = 0.251 ± 0.007 M Q . 

The mass function is often used to determine the mass 
of the compact object, and thus requires estimates of the in- 
clination and companion mass. These parameters are how- 
ever difficult to measure. Several different methods have 
been u sed to estimate the mass of the donor. [Gies fc Boltonl 
(1986) used the spectroscopic orbit, light curve, photo- 
spheric line broadening and an assumed degree of Roche 
lobe filling to find a mass of Mi = 33 ± 10 M©, giving M co = 
16 ±5 M©. From spe ctroscopic analysis of the line spectrum, 
iHerrero et al.l (| 19951 ) derived mass estimates of Mi ~ 18 Mg 
and M co ~ 10 M©. Observa tional and evolutionary con- 
straints led IZiolkow ski (2005) to significantly higher mass 
estimates: M 2 = 40 ± 5 M Q and M co = 20 ± 5 M Q . This 
wide range of masses illustrates the difficulties inherent in 
these indirect measurements of the black hole mass. 

Recently, direct determination of the mass of the com- 
pact object has been suggested from measurements of the 
X-ray radiation. For example, IShaposhnikov fc Titarchukl 
(2007) study the correlation between low- frequency quasi- 
periodic oscillations (QPOs) and spectral index and deter- 
mine the mass of the compact object to 8.7 ±0.8 Mq, which 
is a much stronger estimate than achieved by the indirect 
measurements above. 

In this paper, we will use the value of M co = 8—10 Mq. 
This range is in the lower end of the values found by studies 
of the donor star, but nevertheless compatible with both 
these measurements and the more direct mass estimates of 
the compact object. For the inclination, values from ~ 20° 
to ~ 70° have been suggested, with most falling in the lower 
end of the range. W e will here use i = 28°-38°, following 
IGies fc Boltonl 1 1986ft . These values then give us the current 
system parameters, presented in Table [T] 



2 THE ORBITAL HISTORY OF CYG X-l 

We now turn to the parameters of the system at the time 
of collapse. Our starting point is the currently observed sys- 
tem, summarized in Table [T] The measured value of the 
black hole spin is markedly different from the value inferred 
assuming tidal locking of the binary with its current orbital 
parameters. Assuming a symmetric stellar collapse event, 
this suggests that either the black hole has been spun up 
since its formation, or alternatively that the binary system 
was much tighter in the past, providing a higher angular 
momentum budget for the spinning stellar core at the time 
of collapse. 

To determine the state of the system at the time of 
core collapse we evolve Cyg X-l backwards from its current 
state. In our backwards extrapolation we assume that the 
current mass of the black hole (M co = 9 Mg) is equal to 
the mass of the stellar core which collapsed to create it. 
This is consistent with the picture of very low mass loss, 
less than ~ 1 Mg , implied b y the low velocity of the system 
(|Mirabel fc Rodriguej [20031 ). 

The moment of inertia of the helium core is given by 
/ = fc 2 M co i?co- Before the collapse of the helium core into a 
black hole, the spin parameter a, can be expressed as 



GAfi 



GMco 



(2) 



where Q, co is the rotational frequency and R co the ra- 
dius of the helium core. Assuming a polytropic equation 
of state for the core with index n = 3, gives k — 0.275 
|Lai et al J 11993ft . The radius of a helium co re is given by 
7?co = O.22(M CO /M ) - 6 R© dLee et all 120021) ; a 9 M Q core 
has the radius 0.82 Rq. Knowing the spin and mass of the 
core at the time of collapse thus enables us to determine the 
rotational frequency to be 7.9 x 10 Hz. Under the assump- 
tion of tidal locking this is the same as the orbital frequency 
of the system, corresponding to a period of ~ 0.9 days. We 
will now study the evolution of the system to determine the 
other parameters needed for such an orbit. 



2.1 Evolution by stellar wind 

The system today is wind-fed, and as a first scenario we 
consider the evolution of the binary after the creation of 
the black hole to be driven by mass loss from the donor 
star. This mass loss acts to increase the orbital separation 
while simultaneously decreasing the total mass of the system 
(assuming the accreted mass can be neglected). This implies 
that the binary was tighter in the past than is measured in 
the present. We can estimate this by evolving the binary in 
time. Assuming the fraction of the mass loss accreted by the 
compact object to be much less than unity and that the lost 
mass leaves with specific angular momentum equal to that 
of the do nor, the produc t of separation d and total mass are 
constant (|VerbuntJ| 19931 ) giving 



, Mco + M 2 , , 
d - M co + Mi d ° 



(3) 



where d is the separation and subscript denotes present 
values. The orbital frequency f2 or b of the system (assuming 
Keplerian orbits) is given by 
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< G{M co + M 2 
d 3 



G(M C , 



dg(M co + M 2 , ) 



(4) 



Under the assumption of tidal locking, the orbital fre- 
quency of the core is the same as that of the system. Com- 
bining Eqs. [4] and [2] we now get 



cR 2 



G°-5M co di- 5 (M co + M 2 



-k 2 (M co + M 2 ) 2 



(5) 



Maintaining the assumption that mass accretion onto 
the black hole has not appreciably altered the mass or spin, 
the spin parameter measured today is that of the system 
when the black hole was formed. We can thus use Eq. [5] to 
determine the mass of the donor star at this time. The lower 
panel of Fig. [T] shows the orbital separation as a function of 
donor mass for our range of inclinations. The long-dashed 
line is the s eparation requir e d to give the observed value of 
a, found bv lAxelsson et all (|2005l ). 

To explain the current spin parameter, the mass of the 
donor must therefore have been 46 M© at the time the black 
hole formed, assuming M co = 9 Mq (solid line). Equation[3] 
then gives the orbital separation at this time, 15 Rq. The 
dashed lines show the evolutions allowed when including the 
full range of values in Table [T] Current donor star masses 
of 7 Mq and 21 Mq may at black hole formation have been 
29 Mq and 63 Mq respectively, with the orbital separation 
in these extreme cases being 11 Rq and 19 Rq. 

In order for tidal locking to be viable as the source 
of spin, the tidal synchronisation timescale must be rela- 
tively short. In particular, the tidal locking timescale must 
be comparable to the widening of the orbit as the companion 
loses mass. Using this criterion, we estimate the separation 
at which unlocking occurs by setting the tidal synchronisa- 
tion timescale equal to M/M. We follow the c alculations of 
the sy nchronisation timescale in section 2.3 of iHurlev et all 
(2002), assuming a radiative envelope for the He core. As- 
suming a mass loss rate of ~ 10 -6 Mq yr _1 (see further in 
Sect. 12.41 below), the separation at which the core unlocks 
(and hence the maximum separation at which tidal locking 
can be effective) is ~ 13 Rq. Given the uncertainties in es- 
timating the tidal interactions, this is entirely compatible 
with the orbital separation required by the spin parameter. 

To judge the feasibility of these results, we study the 
ratio between the donor Roche lobe and its radius (upper 
panel of Fig. [TJ , assuming the donor is on the main sequence. 
As the donor mass increases, so does its radius, while the 
orbital separation and Roche lobe radius decrease. For the 
evolution described above, we see that the donor would fill 
its Roche lobe before the required orbital separation and 
mass are reached. In this case, heavy mass transfer would 
occur from the donor to the compact object, and the re- 
sult would likely have been a merger rather than the binary 
required to produce the system we observe today. 

These calculations were performed assuming a poly- 
tropic index of 3, giving a value of k — 0.275. While such 
a polytrope accurately describes massive white dwarfs, it is 
not clear that it correctly matches the stellar core. As this 
most likely contains a central, denser Fe core it may be more 
centrally concentrated than an n = 3 polytrope. Thus, we 
expect that the adopted value of k is on the high side - the 
real value may well be lower. Eqs [3] and [5] show that a lower 
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Figure 1. Orbital separation (lower panel) and donor radius over 
Roche lobe (upper panel) as a function of donor mass when the 
black hole was formed. The long-dashed line in the lower panel in- 
dicates the separation requ ired to give the curren t spin parameter 
of the black hole found bv lAxelsson et all i2005l ). 



value of k leads to an even tighter system being required at 
black hole formation, aggravating the difficulty. 

Clearly, merely adding mass to the donor cannot pro- 
duce the system required at tidal locking. We must therefore 
look more closely at the evolution prior to the formation of 
the black hole. 



2.2 Rotational mixing 

In tight systems, th e rapid rotation of th e stars can lead 
to rotational mixing. Ide Mink et all (|2009l ) recently showed 
that this can lead to chemicaly homogeneous evolution of 
the stars, and the stellar radius will then be much smaller 
than predicted by normal stellar evolution. While they have 
yet to model a syst em such as the ones described above, 
Ide Mink et "all (2009) point to a new channel of binary evolu- 
tion where Roche lobe overflow and subsequent mass trans- 
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fer is significanty delayed or avoided althogeth er due to rota- 
tional mixing (Case M in lde Mink eFaLh oOO). Interestingly, 
tight binary systems with massive stars have been discov- 
ered, e.g. CQ Cep. In this system the stellar masses are 
24Mq and 30Mq, and the orbital period is 1.6d. One might 
therefore speculate whether Cygnus X-l is the result of such 
an evolutionary channel. 

In the previous section we derived the required param- 
eters when the black hole was formed. At this time, the sys- 
tem consisted of a helium core of ~ 9 Mq orbiting a stellar 
companion of ~ 46 Mq. The orbital separation was 15 Rq, 
and the orbital period 0.95 days. The separation can be com- 
pared to the Roche lobe of the donor, which was 8 Rq. This 
is somewhat smaller than the predicted radius of a zero-age 
main sequence (ZAMS) 46 Mq star (which is ~ 10 Rq), and 
comparable to the radius of a star in the lower end of our 
derived range. 

From these results it is clear that even at ZAMS, the 
donor would have been too large. While rotational mixing 
may keep a star from expanding as it evolves, it does not 
reduce the size at ZAMS and cannot be invoked as a means 
of enabling the tidal locking required. 



more realistic value is e = 0.1, giving M cap ture ~ 10 8 M Q 
yr" 1 . The age of the system is difficult to determine, but 
iMirabel fc Rodriguesl l|2003l ) estimated it to be ~ 7 x 10 6 
years by comparing the relative velocity and distance be- 
tween Cyg X-l and its assumed formation site in the Cygnus 
OB3 association of massive stars. The black hole mass could 
thus only have increased by ~ 0.1 Mq in the time since for- 
mation. 

In order to significantly change the spin, the black hole 
must have accreted > 1 Mq. For a 10 Mq black hole, 
the Eddington accretion rate is ~ 10~ 7 Mq yr" 1 . Thus, 
the black hole would need to accrete at the Eddington rate 
for ~ 10 7 years to significantly affect the spin. This is not 
a likely scenario. In addition, the current accretion rate is 
merely a few per cent of the Eddington rate. Therefore, the 
accreted mass can be neglected to first order. The low mass 
accretion onto the black hole also means that the spin pa- 
rameter cannot have been changed by mass transfer from 
the companion star. Ao* < 0.02 for an accreted mass of 
0.1 Mq, hence the value measured now is the same as when 
the black hole was formed. 



2.3 Different evolutionary scenarios 

In order to allow for a wider orbit, and slightly less mas- 
sive companion, one may investigate the results if a slightly 
more massive He core at tidal locking is assumed. The sepa- 
ration is a slow function of mass, so the orbit does not widen 
much. However, the Roche lobe radius is a slightly stronger 
function of mass so the Roche lobe radius of the companion 
actually shrinks slightly even though the orbit widens. It is 
thereby clear that a more massive He core will not allow for 
tidal locking at the required frequency. 

Another alternative to consider is that the system forms 
with two main sequence stars close in mass. The system un- 
dergoes common envelope evolution, leading to two close He 
stars at which point tidal locking occurs. Winds then widen 
the binary to present parameters. The scenario is attractive 
as the secondary does show chemical peculiarities and is en- 
riched in He. However, with a mass ratio close to unity the 
evolution of the two stars would be very similar, and by the 
time the primary becomes a black hole the secondary would 
have lost too much mass to represent the current system. 

2.4 Possible spin-up? 

So far we have neglected any effect of accreted mass and 
angular momentum on the evolution of the system, and 
assumed that the spin is solely determined by the system 
parameters at tidal locking. We can however estimate the 
possible change in spin contributed by accretion onto the 
black hole. The cu rrent mass loss ra te from the donor has 
been esti mated bv Ides et al.1 (|2003h to M wind = 3 x 10" 6 
M Q yr' MVrtilek et all (|2007T ) find a higher mass loss rate, 
M wind = 5 x 10" 6 M yr" 1 . The mass accretion rate is 
clearly variable; however, following IVrtilek et all l|2007l ). we 
can estimate it for a relatively bright state (Lx = 10 38 ergs 
s _1 ). In this case Mcapture = Lx/ec 2 , where e = 0.42 for a 
maximally rotating black hole l|Shapiro fc Teukolskvl[l983T l. 
This yields M capture ~5x 10" 9 M yr" 1 . For Cyg X-l, a 



Effects such as beaming may lead us to underestimate 
the luminosity, and thereby the accretion rate. However, 
studies of X-ray binaries show that the transition between 
the so-called hard and soft states typical ly occur around a 
few per cent of the Eddington luminosity l|Done et al.ll2007r i . 
Spectral studies show that Cyg X-l is mostly observed in the 
hard state; thus, the accretion rate cannot be very high even 
if beaming leads us to underestimate it. Another possibility 
which may lead us to underestimate the accretion rate is 
asymmetric accretion. In this case, the radiative efficiency 
can be much lower than during spherical accretion. While 
such a scenario is unlikely at present, we cannot rule out 
that the system has undergone a period of asymmetric ac- 
cretion in the past, and thereby accreted more mass than the 
estimate above. However, the accretion rate is still limited 
by the capture rate of the stellar wind from the companion, 
and we find it unlikely that asymmetric accretion could have 
increased the accreted mass to the level required. One pos- 
sibility to do so is to invoke a period of super-Eddington ac- 
cretion after the black hole has formed. Such a period would 
have to be very short, and the accretion rate therefore ex- 
tremely high. While this scenario cannot be ruled out, the 
feasibility of super-Eddington accretion is still under debate. 

It is possible that the system underwent a period of 
mass transfer before the creation of the black hole, and this 
may have spun up the stellar core. We will therefore consider 
the case of mass accretion onto the helium core. Although 
mass transfer may have occurred also at earlier stages in 
the binary evolution, accretion at this stage will primarily 
spin up the envelope of the star, and it is unclear how much 
angular momentum is actually gained by the core. 

The maximum time during which the He core can ac- 
crete matter is determined by its evolution, and is of the or- 
der of 10 6 years. In order to increase the rotational frequency 
sufficiently to change the spin parameter by Aa, =0.1, the 
core would have to accrete more than 1 Mq (assuming a 
radius of R co = 0.8/Zq). 
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System 


Mbh 




-^orb 


Reference 




(M q ) 


(Mq) 






Cyg X-l 


9 


14 


5.6d 


see Table [T] 


Cyg X-3 


~ 30 


< 60 


4.8h 


1 


LMC X-l 


11 


32 


3.9d 


2 


LMC X-3 


12 


40 


1.7d 


3 


M33 X-7 


16 


70 


3.5d 


1 


IC 10 X-l 


28 


35 


1.45d 


5 


NGC 300 X-l 


17 


21 


1.35d 


6 



Table 2. Parameters of known HMBHB systems. The na- 
ture of the compact object in Cyg X-3 is still unclear, but 
the curre ntly favored scenario i s tha t of a black hole . Refer - 
enc es: (11 iHia marsdotter et alJ d200E ), (2) lOrosz et all ll2009h. 
(3llYao et all <2005h, (4)lQrosz et al.l J2007h . (5) iPrestwich et al.l 



20071 ). (6) ICrowther et al] |201Ch 



3 OTHER SYSTEMS 

In this context it is also interersting to note that many high 
mass black hole binaries discovered so far have very short 
orbital periods. For comparison, we summarize their param- 
eters in Table [2] 

From these values it is clear that Cyg X-l is not an 
unusually tight system, nor is the spin parameter partic- 
ularly high. The only other high-mass black hole binaries 
with spin estimates are LMC X-l and M33 X-7, but in these 
cases the spin p arameter is estim ated from spect ral fitting 
to be a. ~ 0.9 (|Gou et all 120091 ) and a, ~ 0.8 \hiu et al 



2008), respectively, which is much higher than in Cyg X 
1. Such a high spin is highly unlikely to arise due to tidal 
locking alone. Rather, it would appear that explaining the 
spin parameter in many of the other systems is an even 
greater challenge! The fact that all measured spins are quite 
high reduces the probability for a very peculiar evolution- 
ary channel. In our view, it strengthens the case for a high 
spin being created during the stellar collapse, e.g. through 
asymmetric accretion. For example, it has been shown that 
a spherical accretion shock instability (SASI) can form in 
supernov ae, leading to accretion of sig nificant angular mo- 
mentum (jBlondin fe Mezzacappal l2006). 



4 SUMMARY AND CONCLUSIONS 

We have investigated the possibility for the spin of Cygnus 
X-l to result from the collapse of a stellar core, rapidly spin- 
ning through tidal locking. We find that such a scenario is 
not compatible with standard stellar evolution models; the 
required orbit is so close that the companion star would not 
fit at ZAMS. Rotational mixing or alternative evolutionary 
paths do not offer a solution, nor is the black hole likely to 
have been spun up due to accretion. Our conclusion is thus 
that the spin originates in the black hole formation event. 
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